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Massless QED3 provides a theoretical laboratory for studying infrared (IR) divergences, 
and for investigating dynamical breaking of chiral symmetry [H[2i[3jlU[5j[6j[71[E]. It 
has also emerged as a leading paradigm for understanding the pseudogap phase of cuprate 
superconductors P, [H)]. Lattice simulations of this theory also form an active area of 
research [TTJ [12]. Recently, we have obtained a rather detailed understanding of the IR 
behaviour of this theory [T3l [HI US]- For a particular value of the gauge parameter in a 
non-local gauge [HI (TTJ [HI [HI [20], the IR behaviour is of a conformal field theory with 
canonical dimension for the fermion but a scale dimension one for the photon [13]. This 
value of the gauge parameter is fixed at each order in 1/N (N being the number of fermion 
flavours). For other values of the gauge parameter in this non-local gauge, the IR behaviour 
is again a conformal field theory, but with an anomalous dimension of the fermion which 
is linear in the gauge parameter, and can take arbitrary values. In the usual local gauge, 
the fermion correlation functions are exponentially damped in the IR, as if a mass gap 
is present. However, this is only a gauge artifact and illustrates the pitfalls involved in 
summing the severe IR divergences of the perturbation theory [H] . 

This emphasizes the need to study the IR behaviour of gauge-invariant correlation func- 
tions, as they correspond to physical observables. The simplest such objects are the Green 
functions involving only photons. They have a power law behaviour in the IR, as per a scale 
dimension two (in contrast to the canonical dimension 3/2) for the field strength F^ v (x). 
Another class of objects, which can be related to experiments which probe single fermions, 
are the correlations of gauge-invariant dressed fermions. We have argued in Ref. [15] that 
only the isotropic (in Euclidean space-time) dressing is relevant, and with this dressing, 
the fermion correlations have again a power law behaviour but with a negative anomalous 
dimension. 

There is another class of correlations that are directly relevant to the experiments. These 
are the response functions of QED 3 [21 J corresponding to correlations of gauge- invariant 
composite operators of the form ^}(x)Ti[}(x), where the matrix V is any element of the 
algebra of gamma-matrices. We address their IR behaviour in this paper. We will first 
perform an 0(1/N) calculation which will give an IR logarithm, heralding an anomalous 
dimension for the composite operator. We will then explicitly demonstrate how the ladder 
diagrams lead to an exponentiation of this IR logarithm^] 

A crucial issue in QED3 is dynamical breaking of chiral symmetry and the critical 
number of flavours below which this takes place [U El El [71 [8] • It is always very difficult 
to calculate the condensate < tpip >■ There are different claims for the critical number of 
flavours N c , ranging from N c = 1 to N c = 40 In the course of our analysis, we calculate the 
fermion-antifermion scattering amplitude^] in the scalar channel for vanishing total energy- 
momentum. We demonstrate instability in the vacuum in which this scattering amplitude is 
calculated, for number of flavours less than the critical value N c = 128/37T 2 . Our calculation 
is important for the following reasons. Firstly, we need only ladder diagrams and not self- 
consistent solution of the gap equation. Secondly, we relate the instability to a robust 
physical mechanism: an attractive inverse-square potential (with momentum variables in 

1 The anomalous dimensions for the cases r = 1 and r = 75 in the ladder approximation were also 
addressed in Ref. [Hj. The result is discussed Ref. [25] . 
2 For a review, see, for example, Ref. [21] . 

3 The fermion-antifermion scattering amplitude was also considered in Ref. [25 . 
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the place of coordinate variables) arising from the masslessness of the fermions and the 
feature that the small-momentum behaviour of the photon propagator is inversely linear. 

Let us consider (the Fourier transform of) the correlation function < ip(x)rip(x)ifj a (y)ijj/3(z) > 
(see Fig. [QA). This correlation is gauge- dependent, but from it, one can extract the (gauge- 
invariant) anomalous dimension of the composite operator ipTip, as follows. 

The 0(1/N) contribution to this correlation function, as shown in Fig. [TJ3, is given by 

d?k 1 1 5 ap — £,k a kp/ k 2 



(2vr) 3 la Hi + f f+fi k 2 + pk 

with /i = Ne 2 /8 (for A" four-component spinors). Here we have used the photon propagator 
[T3] in the 1/N expansion (in which there is a resummation of chains of one- loop vacuum 
polarization diagrams on every photon propagator) with a non-local gauge-fixing term0 
This non-local gauge ensures that the IR behaviour of the propagator is inversely linear in 
momentum for arbitrary value of the gauge parameter £. The value of the gauge parameter 
can then be fixed to each order in 1/N such that the IR logarithms are absent in the 
fermion self-energy and in other Green functions of elementary fields [13] . This choice of £ 
will simplify the ladder diagrams, to be considered later in this paper. 

We will now follow the calculation of Appendix A of Ref. [13]. Choose Z M = pLp and 
% — pQn- Also let kp = pKp. Then the integral equals 

c2 f d3K ,. 1 r 1 ,. **, ~ ^k p /k 2 

For p —>■ 0, with and of 0(1), this is IR finite but logarithmically ultraviolet (UV) 
divergent. Letting the divergent part be Clnp, the coefficient C can be obtained by the 
action of [p(d/dp)] p= o on the last expression (it is convenient to replace K by k/p before 
setting p = 0): 

2 f d 3 k ft ft dap — £k a k p /k 2 

C = ~ e Jj2^^¥ 2T ¥ 2lp (k + p) 2 ■ (3) 



4 The gauge- fixing term is given by Eq. (25) of Ref. [M]. The gauge parameter a of Ref. [14] is related 
to the gauge parameter £ of the present paper and of Ref. [13j by a — 1 — £. 
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Figure 2: Ladder diagrams contributing to the correlation function. 



In the numerator, the part proportional to £ equals — £,k 2 T. Writing the other part as 
kixkvlolp^lvlo-, we replace k^k v by |fc 2 <5 M ;,. Now consider two types of vertices [21] : 

lal^l^l* = kT (4) 

with k = 9 if T commutes or anticommutes with all the three gamma-matrices, and k — 1 
if T anticommutes with one or two of the gamma- matrices and commutes with the rest. 
Thus the numerator equals fc 2 (|ft — £)r, and on doing the angular integration, 

Now, if we choose the value £ = 1/3 for which the 0(1/N) logarithmic divergence in the 
fermion self-energy is absent [131, the coefficient of T in C gives the anomalous dimension 
of the composite operator t/Tt/'Jj Thus for the vertices with k — 9, the composite operator 
(for example, ipi/j) has an anomalous dimension 

^ = -3^iV' (6) 

while for the vertices with k — 1, the composite operator does not have an anomalous 
dimension. Our result agrees with that of Ref. [21]. The apparent discrepancy of a factor of 
—2 between Eq. (jH]) and the result of Ref. [2T] is merely due to difference in the definition 
of anomalous dimension. In particular, the authors of Ref. [2T] consider the dependence of 
< ip(x)Tip(x)ip(y)Yip(y) > on the external momentum, and this correlation has two factors 
of the wave-function renormalization Z^ r ^ instead of one factor as in our correlation. It 
is to be noted that we have obtained the anomalous dimension by performing a one-loop 
calculation, instead of the two- loop calculation of Ref. [21j . 

The anomalous dimension was heralded by the appearance of an IR (p — > 0) logarithm in 
external momenta. In order to address the exponentiation of this logarithm, we consider the 
contribution of the ladder diagrams (see Fig. [2]). We choose to work with the special gauge 
parameter £o of our non-local gauge for which the logarithmic terms are absent at all orders 
in the fermion self-energy and in other Green functions of the elementary fields ip, ip and 
A^. Then the full fermion propagator has the behaviour of the free propagator 1/ft, and the 



5 This follows from the application of the Callan-Symanzik equation to the correlation function under 
consideration; see Ref. 1261 for a similar situation. 
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Figure 3: The n-loop ladder diagram: the spinor algebra is drastically simplified for q = 0. 



full vertex cy^, when all the momenta are small. Therefore, keeping only the bare fermion 
propagators and the bare vertices in Fig. [H is justified. Moreover, we will demonstrate the 
exponentiation of only the lowest order contribution to the anomalous dimension. So we 
will not consider the diagrams with crossed photons, which will contribute at higher orders 
in 1/JV. 

Power counting of the subdiagrams shows that there are no UV divergences. We are in- 
terested in the IR behaviour. We will also consider the case q = 0. To extract the anomalous 
dimension of the composite operator, we need the correlation function at non-exceptional 
and Euclidean external momenta. There is a possibility of additional IR divergences at 
exceptional values. However, we consider this case as there is an enormous simplification in 
the spinor algebra, so that the power law in the IR comes out in a simple and straightforward 
way. 

We take T to commute or anticommute with all the gamma-matrices, which is the case 
for which we obtained a log divergence in our one-loop calculation. Now from the Zi-loop 
in Fig. [3j we get the spinor factors 

e*iMn»(&v» - ~ h)&i - k)u/{h - hf) = e 2 (3 - . (7) 

As four gamma-matrices are involved, this contribution is insensitive to our sign convention 
and the choice of the action. Similar simplification occurs for each of the loops of Fig. [3j 
Thus, all spinor dependence is absent, as if the fermion is replaced with a scalar. We get a 
contribution 

d% r d% r dH n 1 1 1 



(27r)»i? J (2^)3(1 J (2n)Hl ^ - l 2 \ ^ - f 3 | ^|f„ - (| 



from n ladders, where we have kept the leading IR behaviour of the photon propagator, 
and presumed a cut-off in the UV for each integration. Integration over l\ leads to ~ In I2 
for small I2 (this is because I2 serves as the IR cutoff for an otherwise log divergent l\ 
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Figure 4: Integral equation for the fermion-antifermion scattering amplitude in the ladder ap- 
proximation (and at vanishing total energy-momentum). 



integration). This when fed into the I2 integration, gives ~ rf(ln/ 2 )ln/2 with an IR cutoff 
Z3, which is ~ ^(ln/3) 2 . Continuing in this way, the n-loop integration of Fig. [3] gives a 
contribution ~ ^j(lnl) n . A sum over the number of ladders n , as in Fig. [21 then gives a 
power in the fermion momentum /, namely, exp[c(e 2 / /i)(3 — £0) In/] = / c ( e2 /^)( 3 ~£o) ; where c 
is a numerical factor. This argument thus suggests an anomalous dimension for composite 
operators like ipip. 

It should be noted that it is the logarithm from the Zi-loop which eventually gets expo- 
nentiated. The vertex in the Zi-loop is not a ip/iip vertex, and so our proof [TBj of absence 
of anomalous dimension for Green functions involving elementary fields in the gauge £ = £0 
does not apply for this correlation function. 

We now analyse this situation using an integral equation. Let us consider the kernel 
1(1, /') depicted in Fig. HI By joining the V lines, we can recover the correlation function 
considered earlier. We have 

1(1,1') = -J^r + A / dH"^ 1 -^- (9) 

I -I' J I" 2 I - 1" 



where 



e 2 



A = 7 ^ r (3-£o). (10) 

Eq. ([9]) clearly holds for the case when the I' lines are joined to give a ipip vertex (that is, 
T = 1), which carries zero spin. Thus, Eq. (Q is for the scalar channel of the fermion- 
antifermion scattering amplitude. We may convert this integral equation into a differential 
equation by operating on it with V 2 : 

(vf + ^)/(Mo = -w 3 )(r-o. (ii) 

This means that 1(1, V) is a propagator (Green function) for the three-dimensional Schrodinger 
equation (with momentum variables in the place of coordinate variables) with the potential 
V(l) = —4n\/l 2 . As A > (£ being 1/3) for us, this is an attractive inverse square poten- 
tial. This is a well-studied problem in quantum mechanics [27], and has unusual properties 
due to exact scale invariance. We will return to this connection a little later. 
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The Green function can be determined using standard techniques [28J. The "radial" 
part fi(l,l') for the channel with angular momentum L satisfies 

Id 2 L{L + l)-47i\ f 4vr 

ldP m P h = ~¥ S{1 (12) 

Note that the the A-containing term (with A > 0) provides a centripetal attraction in con- 
trast to the centrifugal repulsion provided by the angular momentum. The radial equation 
(for I ^ I') has the two power law solutions: 



\ 




4ttA. (13) 



If 47rA < 1/4 (which is true when the number of flavours N > N c = (16/7r 2 )(3 — £o))> a is 
real for each of L — 0, 1, 2, • • •. Now the boundary conditions on fi(l, /') follows from the 
IR and UV finiteness of of the integral on the right-hand side of Eq. (Q: lim^o lfi(l, I') = 
and lirn^oe f L (l, I') = 0. Also, f L (l,l') = f L (l'J)- Therefore, we have f L (l,l') = 
C/ < 1 ^ 2+a / > 1 ^ 2_a , where Z< (/>) is the smaller (larger) of I and I', and C is determined 
from the effect of the delta function in Eq. (|T2|) . Thus, the solution for the Green function 
is obtained as 

nu 7 ) = E E t= U9 A , /, , r nM^yLM^), (14) 

L=0Af=-L ^/(L + |) 2 - 4ttA VM> \«>/ 

where (0, 0) and (#', 0') are the angular variables for Z and I' respectively. 

For the quantum mechanical problem of a particle in an attractive inverse-square po- 
tential, a coupling constant 47rA exceeding the critical value 47rA c = 1/4 leads to a sin- 
gular situation: the particle falls to the centre [27] . In our problem, this corresponds to 
N < N c = 128/37T 2 m 4.3 for £ = 1/3. This is same as the value for the critical number 
of flavours for dynamical chiral symmetry breaking, as obtained from an analysis of gap 
equation [T6], I2"U] . 

Let the eigenfunctions and eigenvalues of the Schrodinger operator —V 2 — AnX/l 2 be 
denoted by Xni}) and e n respectively. Then, from Eq. (iTTj) . 

J(/;/ 7 ) = 4 7 rE Xn(/)X " (//) ■ (15) 

n e n 

There are no bound states (that is, e n < 0) for A < A c [27l[29j, but bound states appear at 
the critical coupling. 1(1, 1') is the Bethe-Salpeter amplitude in the channel with total spin 
zero and total energy-momentum zero. This is possibly the simplest object which brings out 
instability in the vacuum. For total energy-momentum zero, only the mass-squared terms 
in the poles in the amplitude survive. So having some e n < corresponds to a tachyon. 
This means an instability in the vacuum in which we have computed 1(1,1'). The system 
cures this by forming a condensate. 
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Let us now join the V lines in Fig. H] to obtain the desired correlation function. Thus, 
we are to evaluate 



(16) 

in the limit I — > 0. Because of the integral over the angular variables for /', only the L = 
term is picked out from 1(1,1'). Then 

m = *3 - 6Q d[ ,l (LA ^ 



4vrWl/4 -4ttA 



(The full l/(k 2 + fik) photon propagator provides an effective UV cutoff of 0(p).) On 
splitting the range of integration into (0,/) and (/,//), and putting I' = Ix, we get the two 
integrals: 

J' dx x -V2+VV4-4,rA + fa (lg) 

While there is no divergence from the x = end, there is a divergence at x = n/l for Z — > 0. 



On expanding J 1/4 — 47rA in powers of 1/iV, it is found to equal 1/2 + 77 at 0(1/N), where 
?7 is given in Eq. ([6]). Also, the prefactor in the integral in Eq. (1171) equals — rj, to the lowest 
order in 1/N. Thus, 

tufi ( I \ v 

J(l) = -7] dK- 1 -'- - (19) 



for I — > 0. (Since we have a divergence for I —>■ 0, the contribution of unity from the first 
diagram (the free-theory diagram) in Fig. [2] can be neglected.) We have thus demonstrated 
that the correlation function has a power law behaviour with the exponent 

Let us now briefly consider the case q 7^ 0. Let the incoming fermion-antifermion pair 
in the scattering amplitude have the momenta I ± q/2, and the outgoing pair I' ± q/2. Let 
us define L = l/q and L' = I'/q, and scale the loop momentum: L" = l"/q. Also define 
J (I, V , q) = ql(l, I', q). Then J (I, I', q) satisfies almost the same integral equation as Eq. (Q 
for 1(1, 1') with I — > L, V — > L' , I" —>■ L"\ only the factor l/l" 2 in the last term gets replaced 
with 1/ (\L" + q/2\\L" — q/2\), where q is the unit vector along q. (To simplify our analysis, 
we have used 1/k as the fermion propagator.) This then gives the differential equation 




4ttA 



\L + q/2\\L-q/2\ 



^jJ(L,L',q) = -47r<5< 3 > (£-£')■ 



(20) 



So, for large L, we still have an inverse square potential, but near L = ±q/2 the potential 
is Coulomb-like. A potential of this form has also been obtained^ following other routes. 



6 This agrees with the result of Ref. [22] when the appropriate limits are taken. 

7 The authors of Ref. [30] find the same value of N c in QED3 as ours. The two-centre potential appears 
in Section 4 of Ref. ED. 
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In this paper, we demonstrated instability in the scalar channel of the fermion-antifermion 
scattering amplitude in massless QED 3 , for number of flavours less than the critical value 
128/37T 2 for which spontaneous chiral symmetry breaking takes place. This was done using 
only the ladder diagrams, and the instability was linked to the robust physical mechanism 
of an attractive inverse-square potential. We also determined the anomalous dimensions 
of the gauge-invariant composite operators to 0(1/ N) , first as coefficient of IR logarithm, 
and then as exponent of power law. 
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